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Q ■ ABSTRACT 

I We report near-infrared radial velocity measurements of the recently identified donor 

. star in the high mass X-ray binary system OAO 1657-415 obtained in the H band 

■ using ISAAC on the VLT. 

Cross-correlation methods were employed to construct a radial velocity curve with 
a semi-amplitude of 22.1 ± 3.5 km s^^. Combined with other measured parameters of 
this system this provides a dynamically determined neutron star mass of 1.42±0.26 Mq 
, and a mass of 14.3 ± 0.8 Mq for the Ofpe/WN9 highly evolved donor star. 

\^ • OAO 1657-415 is an eclipsing High Mass X-ray binary pulsar with the largest 

I eccentricity and orbital period of any within its class. Of the ten known eclipsing X- 

. ray binary pulsars OAO 1657-415 becomes the ninth with a dynamically determined 

neutron star mass solution and only the second in an eccentric system. Furthermore, 
^ I the donor star in OAO 1657-415 is much more highly evolved than the majority of 

the supergiant donors in other High Mass X-ray binaries (HMXBs), joining a small 
but growing list of HMXBs donors with extensive hydrogen depleted atmospheres. 

Considering the evolutionary development of OAO 1657-415, we have estimated 
the binding energy of the envelope of the mass donor and find that there is insufficient 
energy for the removal of the donor's envelope via spiral-in, ruling out a Common 
\^ . Envelope evolutionary scenario. With its non-zero eccentricity and relatively large 

orbital period the identification of a definitive evolutionary pathway for OAO 1657- 
415 remains problematic, we conclude by proposing two scenarios which may account 
for OAO 1657-415 current orbital configuration. 
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1 INTRODUCTION 
1.1 Neutron star masses 

Despite much theoretical work and the proposal of over 
100 different neutron star (NS) equations of state (EoS) 
l|Kaper at all i2006, ). the precise form of the fundamental 

* Based on observations carried out at the European Southern 
Observatory under programme ID 081.D-0073(B) 
t email: a.mason@opcn.ac.uk 



physical properties of matter under the extreme densities 
and pressures found within NS are still unknown. In this 
regard, observational data can help to eliminate contend- 
ing EoS by rejecting those that place unreasonable con- 
straints on the mass range of observed NS. The recent mea- 
surement of 1.97±0.04 M0 for the mass of the binary mil- 
lisecond pulsar (MSP) PSR J1614-2230 already excludes 
many EoS invoking th e presence of exotic hadronic matter 
iDemorest et al.ll2010l ). while higher NS masses - albeit with 
greater uncertainty - have also been reported (see Section 6). 
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Moreover, one might expect the mass of NSs to show a de- 
pendence on the prior evolutionary history, such that higher 
masses might be expected in binaries where significant accre- 
tion has occurred (e.g. binary MSPs) or where massive pro- 
genitors have y ielded a high pre-SN core mass. With regard 
to the former, I Zhang et al.l (|201ll ) have recently proposed 
that the mean mass of MSPs is greater than that of the to- 
tal NS population, while pTimmes et all (|l996l ) propose a bi- 
modal distribution of NS masses (peaking at 1.27 & 1.76Mq) 
from Type II supernovae (SNe), with a reduced mass for NS 
resulting from Type lb SNe (peaking at 1.32Mq) although, 
as noted by the authors, significant uncertainties in the in- 
put physics remain (e.g. pre SNe mass loss rates). 

Unambiguous NS mass determinations in X-ray binaries 
can only be made for eclipsing binary systems, where the in- 
clination angle is well constrained. In this paper we aim to 
derive masses for both components of the massive eclipsing 
X-ray binary OAO 1657-415, rendered possible by a combi- 
nation of X-ray timing analysis for the NS and radial veloc- 
ity (RV) analysis of the mass donor. The relative scarcity of 
X-ray pulsars, combined with the low geometric probability 
that an eclipse of the binary system will be seen, means that 
only 11 eclipsing X-ray binary pulsars are currently known. 
Up to this tim e 7 of these system s have NS mass determ i- 
nations (e.g. iMason et al.1 (|2010l ). IVal Baker et al.1 (|2005l ). 
iQuaintreU et all l|2003D ~ The NS mass resolution of the re- 
maining systems is thus a priority. 



1.2 The history of OAO 1657-415 

OAO 1657-415 was first de tected over 30 years ago by the 
Copernicus X-ray satellite (|Polidan et al.|[l978l ). Follow up 
observations by HE A 0-1 detected p ulsations with a period 
of 38.22 s (|White fc Pravdd Il979h . whilst later observa- 
tions by the Compton Gamma-Ray Observatory (CGRO) 
BATSE instrument determined that OAO 1657-415 was 
an eclipsing accreting pulsar with an orbital period of ~ 
10.44 days (Chakrabar tv ct al. 1993 ). Along with Vela X-1 
l|OuaintreU et al.ll2003l )^ OAO 1657-415 is one of only two 
eclipsing X-ray binary pulsars to ex hibit a significant eccen- 
tricity, namely e = 0.107 ± 0.001 l|jenke et al.l 120111 ). It is 
further diffe rentiated in poss essing the longest orbital period 
(refined bv lBarnstedt et al.l (,2008 ) as P = 10.448 d) of any 
of these types of eclipsing X-ray binary systems. 

From an initial examination of the orbital parameters of 
the X-ray pulsar it was determined that OAO 1657-415 was 
a high-mass system, indicating the mass of the donor lay in 
the range 14 — 18 M© , with a corresponding radius range of 
25 — 32 R© . Determination of the se stellar parameters led t o 
a spectral classification of B0-6I (IChakrabartv et al.lll993l ). 

The correct identification of the donor in this system 
required the precise location of OAO 1657-415 to be accu- 
rately made. This was achieved by the Chandra X-Ray Ob- 
servatory narrowing the X-ray location error radius down 
to 0.5". Optical imaging of this position did not detect 
any donor candidates down to a magnitude of V>23, im- 
plying many orders of optical extinction for the presumed 
massive mass donor. Near infrared imaging was therefore 
employed, resulting in the identification of a donor located 
within the Chandra error radius. A corresponding IR coun- 
terpart was located in the 2MA5'5' catalogue, namely 2MASS 
J17004888-4139214 with magnitudes J = 14.1, H = 11.7 and 
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Figure 1. Continuum normalised H band spectra of OAO 1657- 
415 (highlighting the Brackctt 11 line) in date order. The orbital 
phase corresponding to each spectrum is indicated on the right 
hand side of the figure. 



Ks = 10.4, an extinction of Ay = 20.4 ± 1.3, and loc ated at 
a distance of 6.4 ±1.5 kpc (jChakrabartv et al.|[20o3 ). 

As a result, near-IR spectroscopy proved mandatory to 
both classify the mass donor in OAO 1657-415 and also 
to construct a n RV curve. Spectr oscopy of the donor ob- 
tained in 2008 (|Mason et al.|[2009l) led to a re-evaluation of 
the spectral classification, with the co rrespondence between 
published NIR BO-6 spectra ((Hanson e t al. 2005j) and that 
of OAO 1657-415 proving poor. Close examination revealed 
that OAO 1657-415 shared a similar spectral morphology 
with that of Ofpe/WNL stars. These are stars in transi- 
tion between the OB main sequence and hydrogen depleted 
Wolf-Rayet stars, whose evolution follows from a wide range 
of progenitor masses. 
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Figure 2. OAO 1657-415 spectral features illustrating the varying wind structure which is particularly noticeable in the P-Cygni line 
profile of the He l 1.700 /im line. 



2 OBSERVATIONS AND DATA REDUCTION 

As the mass donor in OAO 1657-415 is faint (H ~ 11.7) we 
employed the NIR spectrometer ISAAC on the VLT to ob- 
tain high resolution (R ~ 3000) and high S/N spectra in the 
H band. Observations were conducted between 2008 May 
13th and 2008 September 25th in the SW MRes mode with 
a 0.8" slit. Science and RV template exposures centred on 
1.7 fim for a total integration time of 1200s and 128s respec- 
tively were obtained. Spectra were reduced using the ISAAC 
pipeline and wavelength calibration performed utilising OH 
skylines. The resulting data has a count rate below 10 000 
ADU; therefore no correction for non-linearity was required. 

The RV template observed was HD154313, a bright (H 
~ 6.8) standard with a known radial velocity of —29 km s~^ 
and close in spectral type (BOIab) to the previously assume d 
classification of OAO 1657-415 (|Chakrabartv et al.|[l993 ). 
Telluric correction was employed to remove atmospheric fea- 
tures from the target and template spectra. 

Cross-correlation was performed using the standard 
IRAlQ routine fxcor. Unfortunately due to time constraints 
and the end of the VLT observing semester we were unable 
to obtain the full data set of the 20 target spectra we re- 
quested. However, 12 high quality spectra that covered a 
wide range of orbital phase were obtained, sufficient to de- 
termine a dynamical mass solution for OAO 1657-415 (Fig. 
1)- 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



3 DATA ANALYSIS 



After telluric correction the target spectra can be seen to 
exhibit three well defined spectral lines: 1.681 /im Brack- 
ett 11 H I (4-11), 1.700 urn He I and 1.736 /xm Brackett 
10 H I (4-10) (Fig. 2). As noted earlier, OAO 1657-415 has 
been spectrally classified as an Ofpe/WNL star. These tran- 
sitional stars are characterised by exceptionally intense stel- 
lar wi nds with low termi nal velocities and high mass loss 
rates (|Martins et al.|[2007t ). This is demonstrated in the mor- 
phology of the three spectra from the dataset obtained at 
differing orbital phases (Fig. 2), in which P-Cygni profiles 
can be seen in each of the spectral lines, the intensity of 
which clearly varies throughout an orbital cycle. 

Stellar wind contamination leading to varying line- 
profiles make the process of cross-correlation to determine 
radial velocities problematic. After close examination of each 
spectrum within the dataset it was decided to cross-correlate 
around the region surrounding the Brackett 11 1.681 fim 
line. We determined that the measured radial velocities of 
this line would be less affected by the stellar wind, originat- 
ing closer to the photosphere and exhibiting less P-Cygni 
profile variation. Unfortunately for two of the science spec- 
tra (#1 and #12 in Table 1) RV template observations were 
not completed. Additionally 2 science spectra (#3 and #9 
in Table 1) exhibited significantly broadened line-profiles in 
comparison to other spectra in the dataset. To achieve as ac- 
curate an RV measurement as possible it was decided to use 
the science spectrum #2 as a template reference spectrum 
in these four cases. We determined this spectrum possessed 
the most accurate telluric division and the narrowest line 
profile. The resulting radial velocities were then corrected 
to the solar system barycentre and are shown in Table 1. 
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Table 1. The radial velocities of OAO 1657-415 determined from the Brackctt 11 absorption line. 



m number 


Mid-pomt ot Observations (Ui) 


MJD 


Phase 


True anomaly 


Radial velocity / km s ^ 


RV Std 


1 


2008 May 13.340 


54 599.41 


0.249 


0.283 


-45.77 ± 10.3 


Reference 


2 


2008 Jun 09.022 


54 626.22 


0.815 


0.783 


-99.05 ± 10.3 


HD154313 


3 


2008 Jun 18.928 


54 636.00 


0.751 


0.718 


-79.26 ± 10.3 


Reference 


4 


2008 Jun 22.238 


54 639.24 


0.061 


0.075 


-39.93 ± 10.3 


HD154313 


5 


2008 Jun 24.319 


54 641.32 


0.260 


0.293 


-26.02 ± 10.3 


HD154313 


6 


2008 Jun 28.149 


54 645.15 


0.627 


0.605 


-65.87 ± 10.3 


HD154313 


7 


2008 Jul 19.115 


54 666.06 


0.628 


0.606 


-66.72 ± 10.3 


HD154313 


8 


2008 Jul 20.076 


54 667.08 


0.725 


0.694 


-61.13 ± 10.3 


HD154313 


9 


2008 Sept 04.066 


54 713.07 


0.127 


0.153 


-38.67 ± 10.3 


Reference 


10 


2008 Sept 10.054 


54 719.05 


0.700 


0.670 


-79.50 ± 10.3 


HD154313 


11 


2008 Sept 11.013 


54 720.01 


0.792 


0.759 


-81.77 ± 10.3 


HD154313 


12 


2008 Sept 25.110 


54 734.11 


0.141 


0.169 


-45.52 ± 10.3 


Reference 



4 OAO 1657-415 RV CURVE FITTING 

From X-ray data it has been deduced that O AO 1657-415 
has a n elhptical orbit with e — 0.1 07 ± 0.001 ([Jenke et al] 
I2OIII ) and the linear ephemeris of iBarnstedt et al-H 20081 ) 
gives the epoch at which the mean longitude is 90° as 



T90/MJD = 52663.893(10) + 10.44812(13)iV 



(1) 



where the uncertainties refer to the last decimal place quoted 
and A'^ is the cycle number. The mor e recent quadratic 
ephemeris derived by I Jenke et al.l (|201lh concludes that the 
orbital period is currently decreasing with Porb/^orb = 
(-3.34±0.14) X lO"*^ yr"^ but otherwise is very little differ- 
ent from this. Given the relatively low precision of our radial 
velocity measurements, the linear ephemeris is perfectly ad- 
equate for our purposes. During the period in which our 
observations were conducted the accumulated uncertainty 
in phase was calculated as ~ 0.003. 

As OAO 1657-415 has an appreciable eccentricity and 
thus will significantly deviate from sinusoidal motion, we 
cannot fit the radial velocities of the donor in this system 
with a sinusoidal radial velocity curve. For an eccentric sys- 
tem we can construct a sinusoidal radial velocity curve by 
fitting the true anomaly u (the angle between the major axis 
and a line from the star to the focus of the ellipse) to the 
calculated radial velocities of the donor. In order to do this 
we must first deduce the mean anomaly M (the time since 
the last periapsis multiplied by the orbital phase). This in 
turn is related to the eccentric anomaly E (the angle be- 
tween the major axis of the ellipse and the line joining the 
position of the object and the centre of the ellipse) by the 
equation 



2n 

E - esinE ^ M = —{t - To) 



(2) 



The right hand side of the above equation is the orbital 
phase, with To the time of periastron passage. In an eccentric 
orbit. To and Tgo are related by 

P{uj - tt/2) 



To = Tgo + ■ 



2tt 



(3) 



In the case of OA O 1657-415. the l ongitude of periastron, 
w = 91.7° ± 0.5° jjenke et al.ll201ll '). so the offset between 
the two is only about 0.005 in orbital phase. 

This eccentric anomaly equation cannot be solved ana- 
lytically, but solving it numerically we obtained the eccentric 
anomaly for each observation. These eccentric anomaly val- 




0.5 

Itiie anomaly / 2n 

Figure 3. Radial velocity solution for the mass donor in 
OAO 1657—415. The dashed line represe nts a fixed zero p hase co- 
incident with the published ephemeris jBarnstedt et aU 1(2008)). 
The solid line is the best fitting sinusoid with three free parame- 
ters, allowing for a shift in phase. 



ues can then be finally related to the true anomaly by the 
equation 



tan — 
.2 



(i 



l + e 



tan 



(I) 



(4) 



The calculated true anomaly values are shown in Table 1 and 
the measured radial velocities plotted as a function of the 
calculated true anomaly are shown in Figure 3. Over-plotted 
on these data are two sinusoidal fits. The first has just two 
free parameters (RV amplitude and systemic velocity) with 
a period and zero phase fixed by the ephemeris of Barnst- 
edt et al. (2008). The second fit allows the zero phase as a 
third free parameter to allow for any possible accumulated 
uncertainty. However, this is only a marginally better fit and 
the reduction in chi-squared does not justify the inclusion of 
an extra free parameter. We therefore use the first fit which 
yields a radial velocity amplitude of 22.1 ± 3.5 km s"'^ and 
a systemic velocity of 57.2 ± 3.0 km s~^. Error bars on 
each radial velocity value shown in Figure 3 are such that 
the reduced chi-squared of the fit is equal to unity. 
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5 CALCULATION OF SYSTEM PARAMETERS 

Jenke et al (2011) measured the projected semi-major axis 
of the neutron star's orbit from X-ray pulse timing delays as 
Ox sin i = 106.16 ± 0.08 light seconds. From this, the semi- 
amplitude of the neutron star's radial velocity may be cal- 
culated using 



ax sm I ■ 



(5) 



to give ATx = 222.8 ± 0.2 km s"^ 

To determine the masses of the system components we 
must first consider the mass ratio of the system q which 
is equal to the ratio of the semi-amplitudes of the radial 
velocities for each star 

_Mx _Ko 
^ Mo Kx 
where Mo and Mx are the masses of the supergiant star and 
neutron star respectively, with Ko and Kx the correspond- 
ing semi-amplitudes of their radial velocities. Specifically for 
elliptical orbits. 



(6) 



Mo = 



iCx'P(l-e)5 



2tvG sin^ i 
and similarly 



Mx 



Ko'^Pjl-e)^ 
2ttG sin^ i 



1 + 



(7) 



(8) 



where i is the angle between the normal to the plane of the 
orbit and the line-of-sight and P is the orbital period. 

The geometry of the eclipse yields the following relation 
between the angle of inclination of the system, i, the eclipse 
half angle, 0, and various length scales in the system: 



Here Rl is the Roche lobe radius of the supergiant, r is the 
separation between the centres of the two stars and /3 is the 
ratio of the supergiant's radius to that of its Roche lobe, re- 
ferred to as the filling factor. For stars in a circular orbit, Rl 
and /3 are constant, and the separation between the centres 
of the two stars is just the radius of the or bit. An approxima- 
tion o f the Roche lobe radius is given by I Joss fc RappaportI 
(Il984h as 



A-)-_Blogg-(-Clog^5 



where they determine the constants as 
A ^ 0.398 - 0.0260^ + 0.004Q^ 



B ^ -0.264 
C ^ -0.023 



-0.0520^ - 0.015^'' 
0.005^^ 



(10) 

(11) 
(12) 

(13) 



fl is the ratio of the spin period of the supergiant to its 
orbital period. Since the system has not circularized, Q is 
unlikely to be unity. If tides are efficient, then we might 
suppose that the rotation of the donor is synchronized at 
periastron, where the tides are most effective. In this case, 
the star's angular velocity will be larger than the orbital 
velocity at mid-eclipse, and at this point Q ~ 1.2. To account 
for this uncertainty in our calculations we therefore adopt 
Q — 1.2±0.2 as a conservative estimate. 



For systems with a significantly eccentric orbit such as 
OA01657-415, the Roche lobe radius, the separation be- 
tween the centres of mass of the two stars, and the Roche 
lobe filling factor will all vary with orbital phase, although 
Rh/r can be assumed to be constant. Using the above, we 
therefore calculate Rh/r — 0.553 ± 0.025. 

The remaining system parameters were then calculated 
under two assumptions. If the donor star fills its Roche lobe 
at periastron, then the filling factor at mid-eclipse will be 
P ~ 0.9. This allows us to calculate a lower limit to the 
system inclination, i, and upper limits are arrived at for the 
stellar masses. Conversely, if the system is viewed edge-on 
(i — 90°), we calculate a lower limit for the Roche lobe 
filling factor, /3 at mid-eclipse, and lower limits are derived 
for the stellar masses. The limits on the semi-major axis of 
the orbit, a, are simply determined from Kepler's third law 
using the calculated ranges for the masses of the two stars. 
In practice, it turns out that the uncertainties on the system 
parameters overlap significantly for these two limits, so in 
Table 2 we list a single set of values. 

Now, the separation between the centres of the two stars 
at mid-eclipse is given by 

r = (14) 
1 + e cos ui 

where oi is the argument of periastron, measured by Jenke 
et al (2011) to be a; = 91.7° ± 0.5° in this case. This yields 
a value of r = 49.9 ± 1.0 R0. Then the Roche lobe radius, 
_Rl, is this separation multiplied by the factor 0.553 ± 0.025 
calculated earlier, which gives the value listed in Table 2. 
Finally, the radius of the donor star, _Ro is just /3-Rl, and 
this value is also shown in Table 2. 

Uncertainties on each of the inferred system parameters 
shown in Table 2 were determined by a Monte-Carlo calcu- 
lation which propagated the uncertainties in each observed 
parameter through one million trials, assuming a Gaussian 
distribution in each individual uncertainty. 



6 DISCUSSION 

6.1 System properties 

The H band spectra for the period May - July, 2008 (Fig. 2) 
demonstrate a high degree of line profile variability (LPV) 
in the Hi transitions, with the Hi (4-10) line observed to 
transition from absorption to a P Cygni profile. He I 1.70 
^m is also present and displays a variable P Cygni profile, al- 
though He II 1.69/xm is absent. Conversely, the K band spec- 
trum of this object previously published (Fig 2. lMason et al.l 
lioO^), presents the Hell 2.189 fim line in emission. The H 
band spectrum is more reminiscent of P Cygni early B hy- 
pergiants than the Ofpe/WNL classification that had been 
inferred from the K band data alone. Indeed, the H and K 
spectra can only be reconciled if it is assumed that the He ll 
2.189 fim transition is driven into emission by the presence 
of the neutron star. 

Precedents for excess Hell emission in HMXBs in both 
the optical and near-IR are provided by LMC X-3 fc X-4 
jNegueruela fc Coell2002l ) and CI Cam (|Clark et al.lll999l ) 
respectively. Radial velocity shifts in this line in LMC X- 
4 confirm its association with the compact accretor, with 
the observed line strength variability possibly resulting 
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Table 2. System parameters for OAO 1657-415. 



Parameter 



Value 



Ref. 



Observed 

ax sin i / It sec 

P / d 

Tgo / MJD 



Be I deg 

Ko I km s-l 

Assumed 

Inferred 

iCx / km s~^ 

i I deg 

Mx / Mq 
Mo I M0 
a / Rq 
iJL / R© 



[1] 
[1] 
[3] 
[4] 



106.16 ±0.08 [1] 
10.44812 ± 0.00013 [2] 
52663.893 ± 0.001 [2] 
91.7 ±0.5 
0.107 ±0.001 

29.7 ± 1.3 
22.1 ±3.5 

1.2 ± 0.2 

222.8 ±0.2 
0.099 ± 0.016 
0.896 ±0.051 
87.0 ±7.7 
1.42 ±0.26 
14.3 ±0.8 
50.3 ± 1.0 
27.6 ± 1.6 

24.8 ± 1.5 



[1] Jenke et al. 2011; [2] Barnstedt et al. 2008; 
[3] Chakrabarty et al. 1993; [4] this paper 



from changes in accretion rate; similar observations for 
OAO 1657-415 would test the hypothesis that the line is 
likewise associated with the accretor rather than the pri- 
mary. The significant rapid, non-secular LPV in the remain- 
ing H I and He I lines is most easily understood by a depar- 
ture from spherical symmetry by the presence of the neutron 
star. The WNLh star Wdl-44 provides a precedent for this 
conclusion, also demonstrating dramatic LPV in wind dom- 
inated H I and He I P Cygni profiles, although in this case it 
is expe cted that the corn panion is a normal H/He-burning 
object ('Clark et allbOld '). 

Mindful of the above, we attempted to derive the phys- 
ical and wind properties of the donor star in OAO 1657- 
415 using the stellar a tmosphere modelling code CMFGEN 
iHillier fc Miller 19981) and th e meth od as described in Sec- 
tion 4.1 of lCrowther fc Evand (|2009l 'l. utilising the H and K 
band spectra obtained between May - September 2008. Un- 
fortunately, we were unable to obtain a satisfactory formal 
fit to these data, which we attribute to the reasons above, 
namely a lack of contemporaneous data for a star showing 
LPV - possibly as a consequence of an aspherical wind - 
and the presence of an additional source of ionising radia- 
tion within the system yielding discrepant He ll 2.189 /im 
emission. 

While we were unable to obtain a formal fit to the 
model, we were able to determine order of magnitude es- 
timates for the bulk properties of the star that we regard as 
indicative of its true properties: T ~ 20 kK, L ~ lO^'^L©, 
- 30 R©, M = 2 X 10"'^ M© yr"\ ~ 250 km s"\ 
Using our spectra from the H and K bands the model also 
provided us with estimates of the mass fraction of hydrogen 
Xh ~ 5% (H/He = 0.21 by number) in which our primary 
diagnostics were the Brackett series lines (1.680, 1.736 and 
2.166 /im), helium Xhe = 94% by mass from He I lines 
(1.700, 2.058, 2.112, 2.161 fim) and other metals Xz = 1%. 

Comparing the stellar parameters of OAO 1657-415 



with the Ofpe/WNL star AF ^Martins et al."2007'), we find a 
correspondence between temperature and radius, whilst the 
luminosity of OAO 1657-415 [L ~ lO^'^L©) is consistent 
within the uncertainties to that of AF. However, the lower 
than expected luminosity of OAO 1657-415 implies that it 
could be a B supergiant rather than a Ofpe/WNL star. Not- 
ing that there is a large degree of overlap between model 
derived luminosities of Ofpe/WNL and B supergiants, we 
turn to other stellar properties to address this possibility. We 
find the low terminal velocity of OAO 1657-415 combined 
with its mass loss rate and temperature a re incompatible 
with those found for field Blue Supergiants (|Crowther et al.l 
l2006h . B supergiants with a luminosity of L ~ lO^'^L© would 
be expected to have a radius greatly in excess of that found 
from both modelling a nd dynamical calculations of ~ 30 
Mq (jSearle et al.ll2008l ). B supergiants also have a signifi- 
cantly greater hydrog en abundance of H/He = 5 by number 
l|Crowther et al.ll2006l ) than that found during our modelling 
of OAO 1657-415. With the degree to which hydrogen is 
depleted within OAO 1657-415 its radius may seem exces- 
sively large. However, several examples of WNL stars with 
atmospheres with a similar or greater degree of hydrogen de- 
pletion and comparable radii to that of OAO 1657-415 have 
been detailed in the literat ure. In particula r, the emission 
line star HD 326823 tMarcolino et al.l bOOTl ) is reported as 
entering the Wolf-Rayet phase as a WN8 star with a radius 
R = 30 R© and a hydrogen mass fraction of Xh ~ 3%. Simi- 
larly the WN8 star WR123 (|Crowther et al.lll995l ') may have 
a hydrogen mass faction as low as 0.5% yet is reported as 
having a radius R ~ 15 R©. 



6.2 OAO 1657—415 in an evolutionary context 

When these properties are combined with the finding that 
OAO 1657-415 possesses a severely H-depleted atmosphere, 
it is apparent that it is a rather evolved object and is clearly 
more advanced than the typical OB SG mass donors in 
HMXBs. A comparable mass donor in this regard is Wray 
977, the blue hypergiant companion within GX301-2. Prior 
to OAO 1657-415, Wray 977 was the most evolved non-WR 
donor known within a HMXB system. Wray 977 has a sig- 
nificantly higher hy drogen mass fraction of 45%, (Table 8 in 
|KaEeretalJ(|200i)) than that of OA01657-415. 

Thus, while a definitive spectral classification for the 
mass donor in OAO 1657-415 is difficult, even if it has 
yet to formally become a WR star, it is apparent that it 
will do so shortly. As such it will become only the fourth 
known HMXB with such an evolved mass donor. The only 
other example in our Galaxy is Cyg X-3, which has a WNE 
mas s donor orbiting a NS or B H companion in a 4.8 hr or- 
bit (|van Kerkwiik et al.lll996l ). Ivan den Heuvel fc De Loord 
(|i97a r describe an evolutionary scenario for this system, 
whereby the mass donor in the OB SG+accretor progeni- 
tor system expands towards a RSG phase, triggering Com- 
mon Envelope evolution which leads to the formation of the 
WNE star and a dramatic contraction in orbital separation. 

In contrast, while the two known ex tragalactic exam- 
pl es - ICIO X-1 llClark fc Crowthej 12003 ) and NGC300 X- 
1 (|Crowther et al.1 12007^ - also host WN stars, they both 
have significantly longer orbital periods (~ 30 hrs), and ap- 
pear to host more luminous mass donors and more massive 
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accretorfl These properties preclude these binaries from 
evolving via the scheme inferred for Cyg X-3. Instead, it 
is thought that they must have evolved from very compact 
ini tial 0+0 star binarie s via the case M pathway described 
bv lde Mink et al.l (|2009l '). In this, the rapid rotation of both 
(tidally locked) components leads to efficient mixing and 
hence homogeneous chemical evolution, which results in the 
stars remaining within their Roche lobes while evolving to 
WR stars. This evolution prevents binary driven mass loss, 
resulting in a longer final orbital period and crucially a heav- 
ier relativistic companion. 

However, the properties of OAO 1657-415 do not fit 
either of these schemes. The non-zero eccentricity appears 
to preclude post-SN common envelope evolution (as a re- 
sult of circularization processes during spiral-in) and hence 
OAO 1657-415 being the endpoint of the van den Heuvel & 
De Loore (1973) pathway. To address this i n a more quan- 
titativ e manner we adopt the formalism of iDewi &: Taurig 
l|2000l) , whereby the evolution of a putative RSG+NS sys- 
tem through a common envelope (CE) phase may be fol- 
lowed by the paramerisation of the binding energy of the 
envelope of the donor star (A) and the efficiency parameter 
r] for the removal of the envelope during spiral-in. In the 
following discussion we assumed an initial donor mass of 
~ 40 Mq resulting in a ~ 24.5 Mq envelope and ~ 15.5 Mq 
core, representative of the current mass of the Ofpe/WNL 
star. The binding energy of the envelope of the mass donor 
is given by 



£^bind — 



-GM, 



donor -^^env 



o 



(15) 



where Mdo 



Monv and -Rdonor are the masses of the donor 



and envelope of the donor and the donor radius, respec- 
tively. In order to facilitate a successful ejection of the donor 
star envelope we consider a wide orbit prior to the CE 
for two reasons: A wide orbit requires an evolved donor 
star to fill its Roche-lobe and giant stars have relatively 
small binding energies of their envelope. Furthermore, a 
wide pre-CE orbit has the potential to release more orbital 
energy during spiral-in. Even for a rather evolved donor 
(i^donor — 1500 J?0) we find j-Ebindl ^ 10^° e rg giv en that 
for such a star A ^ 0.02 (c.f. jPewi fc Taurig (|200ll 'l'l. How- 
ever, simply determining the available orbital energy; 



A'/ns 



2a 



2ao 



8.2 X lO^'^erg (16) 



where the mass of the neutron star, Mns = 1.5 A/0 and ag 
and a are the initial and final orbital separations. Thus even 
assuming the smallest possible envelope binding energy and 
a 100% ejection efficiency (rj — 1.0) for conversion of orbital 
energy we find that there is insufficient energy available for 
ejection of the donor envelope (by 2 orders of magnitude). 
This is indicated graphically in Fig. 4 which demonstrates 
that for a successful ejection of the envelope yielding the 
current orbital parameters one would require A > 1 (solu- 
tions are only possible above the dashed line), while stellar 
structure calculations for a 40 Mq star yield A = 0.006 — 0.02 
(Dewi & Tauris 2001). 



2 > 15M0 for NGC 300 X-1 fllCrowther et al.l boicl 'l') and > 
23M0 for ICIO X-1 fllSilverman fc Filippenkd l2008t) 'l versus < 
IOM0 for Cyg X-3 ( jHanson et al.lll999h '). 



If the common envelop e evol utionary pathway 
of Ivan den Heuvel fc De Loor j j 19731 ) cannot produce 
OA01657-415 what might? The large degree of H-depletion 
in the mass donor implies a post-rather than pre-RSG 
state; hence it is difficult to interpret the system as evolving 
towards such a common envelope phase. Moreover, case 
M ch emically homogeneous evolution (e.g. Ide Mink et al.l 
Hooi)) also seems inappropriate since it would result in a 
much more massive relativistic companion than we find. 

In the absence of any tailo red numerical model for the 
evolution of OAO 1657-415 fc.f. lWellstein fc Langeil (|l999al ) 
for GX301-2) any scenario is inevitably somewhat specula- 
tive. It might be supposed that the progenitor system con- 
sisted of two O stars in a comparatively short orbit. Mass is 
transferred from the primary to the secondary, which after 
rapidly gaining mass evolves more quickly than the donor 
and hence undergoes SNe first, leaving the primary in a 
H depleted phase. Despite such mass transfer leading to a 
rather wide pre-SN binary configuration, a fortuitous SN 
kick may reduce the separation, with subsequent interaction 
acting to circularise the orbit. While this would leave the pri- 
mary in the requisite H-d epleted state, simulatio ns of conser- 
vative mass transfer (e.g. IWellstein et al.l (|200ll )) leading to 
reverse SNe yield very low masses for the mass donor in com- 
parison to our dynamical mass estim ate. Alternatively, f rom 
our earlier paper on OA01657-415, (jMason et al.ll2009l ) lu- 
minosity calculations a nd comparison to evolutio nary tracks 
for massive stars from iMevnet fc Maedej (|2000|) indicate a 
progenitor mass for the 0A01657-415 donor of ~ 4OA/0. 
This implies OAO 16 57-415 is massive en ough to have under- 
gone an LBV phase (|Langer et al.|[l993 'l. If the evolution of 
the progenitor system had proceeded such that a LBV+NS 
was produced with an orbital period in which both com- 
ponents would interact during the LBV phase, in this case 
the binding energy of the donors envelope would be reduced 
whilst the corresponding efficiency parameter rj for the re- 
moval of the envelope would increase. The changes to these 
parameters might be sufficient to enable the removal of the 
donors envelope by common envelope ev olution. This sce- 
nario has been suggested qualitatively bv iDe Donder et al.l 
(jl997h . The hydrogen depleted HMXB system GX301-2 
which hosts a blue hypergiant donor of 43Mq ± lOAf0 with 
a 1.85 ± O8Af0 NS companion in an eccentric 41d orbit may 
provide an analogue to the pro genitor system to OA01657- 
415. ModeUing of this system bv lwellstein fc Lange"3 l|l999bl ) 
as a 25 -I- 26Af0 binary which undergoes conservative case A 
mass transfer with an initial 3.5 day period was able to re- 
produce the current observed mass of the donor successfully, 
demonstrating that it is possible for a massive hypergiant to 
have a NS companion. 

We are thus currently unable to determine a definitive 
qualitative evolutionary pathway for OAO 1657-415, but 
conclude that it cannot have formed via either of the sce- 
narios proposed for the other three known WR X-ray bina- 
ries - consequently this class of object appear to be rather 
heterogeneous in nature in comparison to other subsets of 
HMXB such as the OB SG and Be star X-ray binaries. 

However we do note that the future evolution of the 
mass donor in OAO 1657-415 will result in an increase in 
wind velocity, while wind driven mass loss will cause the 
separation to increase. The combination of both effects will 
result in a reduction in accretion efficiency and hence X-ray 
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Mass ratio parameter, k (qp/q) 

Figure 4. Plot of the initial/final mass ratio (in this case q is 
defined as A/core/MjMg), fc(= qo/q), versus the change in the or- 
bital period - P/Pq- The dashed line represents the value of the 
latter assuming a maximum orbital separation of ~ 2500 Rq 
(Po = 2250 days), while the shaded region corresponds to the 
range valid for k. The curved lines represent different binding en- 
ergies of the envelope of the donor star (A), assuming g = 10 and 
an ejection efficiency r] = 0.3. See Sect 6.2 



luminosity. Hence one might suppose that even if the evolu- 
tionary pathway that formed O AO 1657-415 is traversed by 
large numbers of stars, comparatively few will be detectable 
as X-ray sources due to the likely rapidity of passage through 
the Ofpe/WNL phase. 



6.3 Future evolution of OAO 1657-415 

The final mass of the Wolf-Rayet star prior to its core- 
collapse is probably somewhere between 5 — 10 Mq 
l|Mevnet fc Maeder, .200^ 1. If the star loses its mass in the 
form of a direct fast wind, (assuming the vast majority of 
the mass is lost from the binary) the system will have an 
orbital period of more than 65 days when the companion is 
reduced to 5 Mq prior to the collapse in less than 1 Myr. If 
the outcome of the core collapse is a neutron star then the 
ejection of the envelope surrounding the iron core is likely 
to disrupt the binary; a binary will be disrupted if it loses 
more than half of its total mass in an instantaneous event 
like a supernova (|Colgatd[T970l ). On the other hand, the 
momentum ki ck expected to be imp arted to the newborn 
neutron star (|Lvne fc LorimeJ Il994l ) can also, depending 
on its direction, have the effect of tightening the binary so 
it survives and remain s bound d espite losing more than 
half of its total mass (|Hillslll983l ). In this case a double 
neutron star system is formed and would most likely be 
wider than the widest of the currently 9 double neutron 
star systems known in our Galaxy, PSR J1811-1736, which 
has an orbital period of 18.8 days. 



7 CONCLUSIONS 

Following on from the recent spectral classification of the 
High Mass eclipsing X-ray binary pulsar OAO 1657-415 
((Mason et al.ll2009l ). we have reported the analysis and re- 
sults from multi-epoch observations of the object using 
the ESO/VLT NIK spectrograph ISAAC. We have con- 
structed an orbital solution from the measured radial veloc- 
ities which, combined with the orbital parameters of this ec- 
centric system, provides a dynamically determined neutron 
star mass of 1.42 ± 0.26 Mq. The donor star in this system 
was previously determined to be of type Ofpe/WN9, that 
has experienced a significant degree of mass loss. The radial 
velocity curve constructed showed that the donor mass lies 
within the range 14.3 ± 0.8 Mq. 

An attempt was made to model the stellar atmosphere 
of the donor; this was found to be problematic due to the 
presence of He ll emission that made a formal fit to the 
spectrum difficult. Modelling of the stellar atmosphere did 
however provide order of magnitude estimates of stellar pa- 
rameters, indicating that the Ofpe/WN9 star in this sys- 
tem, although significantly hydrogen depleted, has retained 
a small fraction of its envelope. This combined with the non- 
zero eccentricity and 10 day orbital period of the system 
is indicative of a different evolutionary pathway than in the 
three other WR-CC systems known. Calculating the binding 
energy of the envelope we found there is insufficient energy 
to eject the envelope. This combined with the orbital pe- 
riod of the system effectively rules out Common Envelope 
evolution as a viable evolutionary scenario. We have pos- 
tulated other evolutionary pathways for OAO 1657-415 in 
which mass transfer occurs in a short binary orbit in which 
the secondary is the first to experience a SNe, or the system 
experiences a 'mild' form of CE evolution during an LBV 
phase. Whilst we are unsure of the exact evolutionary pro- 
cess responsible for the current configuration of OAO 1657- 
415, we are convinced that the evolution of this system forms 
an entirely new evolutionary scenario than the two proposed 
to explain the development of the other three known Wolf- 
Rayet X-binary systems. 
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